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ABSTRACT 


A  rapid,  accurate  method  was  developed  of  computing  propa¬ 
gation  loss  as  a  function  of  range  in  the  ice  covered  Arctic 
Ocean.  Input  parameters  to  the  propagation  model  are  source 
and  detector  depth,  wave  frequency,  ice  roughness,  bottom  topo¬ 
graphy,  and  the  velocity  structure  as  a  function  of  depth  in 
the  icq  water,  and  bottom.  Computation  is  done  by  direct  inte¬ 
gration  of  the  exact  integral  solution  of  the  wave  equation  de¬ 
rived  from  a  harmonic  point  source  located  in  a  multilayered,  in- 
terbedded  liquid-solid  half  space.  Tne  integration  technique  in¬ 
troduced  by  H,  W.  Marsh,  employs  the  Fast  Fourier  Transform  for 
very  rapid  evaluation  of  the  integral  solution.  Computed  propa¬ 
gation  loss  as  a  function  of  range  is  in  good  agreement  with 
field  data. 
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Figure 


CAPTIONS  FOR  FIGURES 


1.  Typical  ray  paths  for  a  source  100  m  deep  in  the  cen¬ 
tral  Arctic  Ocean.  Forty  rays  computed  at  1  degree 
intervals  at  the  source.  Angles  at  the  source  gc 
from  20  degrees  above  to  20  degrees  below  the  horizon¬ 
tal. 


2.  Oscillogram  of  signal  from  9  lb.  charge  of  TNT  fired 
at  a  depth  of  152  m.  Hydrophone  at  a  depth  of  61  m. 
Passband  of  listening  system  10-21  Hz. 


3.  Comparison  of  observed  and  computed  group-velocity  dis¬ 
persion  for  first  three  modes.  Range  578.6  km.  5-lb. 

TNT  charge  at  66  m,  hydrophone  at  46  m. 


4.  Vertical  sound  velocity  profile  for  model  of  Table  I. 


5.  Propagation  loss  as  a  function  of  range  at  10  Hz  for 
model  of  Table  I.  Direct  integration  (FFP) .  Source 
depth  150  m,  hydrophone  depth  100  m. 


6.  Propagation  loss  as  a  function  of  range  at  25  Hz  for 
model  of  Table  I.  Direct  integration  (FFP), Source 
depth  150  m,  hydrophone  depth  100  m. 


7.  Propagation  loss  as  a  function  of  range  at  10  Hz  for 
model  of  Table  I.  Normal-mode  theory.  Source  depth 
150  m,  hydrophone  depth  100  m. 


8.  Propagation  loss  is  a  function  of  range  at  25  Hz  for 
model  of  Table  I.  Normal- mode  theory.  Source  depth 
150  m,  hydrophone  depth  100  m. 


9.  Absolute  value  of  integrand  at  10  Hz  corresponding  to 
Figure  5. 


10.  Absolute  value  of  integrand  at  25  Hz  corresponding  to 
Figure  6. 


11.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  2  at  10  Hz.  Source  depth  50,  hydrophone  depth 
50  m.  Surface  roughnerjs  3  m  RMS. 
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Captions  for  figures  (cont'd) 


Figure  12,  Propagation  loss  as  a  function  of  range  for  model  of 
Table  II  at  15  Hz.  Source  depth  50  m,  hydrophone 
depth  50  m.  Surface  roughness  3  m  RMS. 


Figure  13.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  II  at  20  Hz.  Source  depth  50  m,  hydrophone 
depth  50  m.  Surface  roughness  3  m  RMS. 


Figure  14.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  II  at  30  Hz.  Source  depth  50  m,  hydrophone 
depth  50  m.  Surface  roughness  3  m  RMS. 


Figure  15.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  II  at  60  Hz.  Source  depth  50  m,  hydrophone 
depth  50  m.  Surface  roughness  3  m  RMS. 


Figure  16.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  III  at  10  Hz.  Source  depth  50  m,  hydrophone 
depth  50  m.  Ice  roughness  3  m  RMS. 


Figure  17.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  III  at  15  Hz.  Source  depth  50  m,  hydrophone 
depth  50  m.  Ice  roughness  3  m  RMS. 


Figure  18.  Propagation  less  as  a  function  of  range  for  model  of 
Table  III  at  20  Hz,  Source  depth  50  m,  hydrophone 
depth  50  m.  Ice  roughness  3  m  RMS. 


Figure  19.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  III  at  30  Hz,  Source  -epth  50  m,  hydrophone 
depth  50  m.  ice  -  ughness  3  ir.  RMS. 


Figure  20.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  III  at  60  Hz,  Source  depth  50  m,  hydrophone 
depth  50  m.  Ice  roughness  3  ro  RMS. 


Figure  21.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  IV  at  10  Hz.  Source  depth  50  m,  hydrophone 
depth  50  m.  Surface  roughness  3  m  RM3. 
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CAPTIONS  FOR  FIGURES  (cont'd) 


Figure  22.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  IV  at  15  Hz.  Source  depith  50  m,  hydrophone 
depth  50  m.  Surface  roughness  3  m  RMS. 


Figure  23.  Propagation  loss  as  a  function  or  range  for  model  of 
Table  IV  at  20  Hz.  Source  depth  50  m,  hydrophone 
depth  50  m.  Surface  roughness  3  m  RMS. 


Figure  24.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  IV  at  30  Hz.  Source  depth  50  m,  hydrophone 
depth  50  m.  Surface  roughness  3  m  RMS. 


Figure  25.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  IV  at  60  Hz.  Source  depth  50  m,  hydrophone 
depth  50  m.  Surface  roughness  3  m  RMS. 


Figure  26.  Propagation  loss  as  a  function  of  r.nge  for  model  of 
Table  V  at  10  Hz.  Source  depth  50  m,  hydrophone 
depth  50  m.  Ice  roughness  3  m  RMS.  . 


Figure  27.  Propagation  loss  as  a  function  of  range  for  model  of 
Taule  V  at  15  Hz.  Source  depth  50  m,  hydrophone 
depth  50  m.  Ice  roughness  3  m  RMS. 


Figure  28.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  V  at  20  Hz.  Source  depth  50  m,  hydrophone 
depth  50  m.  Ice  roughness  3  m  RMS. 


Figure  29.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  V  at  30  Hz.  Source  depth  50  m,  hydrophone 
50  m.  Ice  roughness  3  m  RMS. 


Figure  30.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  V  at  60  Hz.  Source  depth  50  m,  hydrophone 
depth  50  m.  Ice  roughness  3  m  RMS. 


Figure  31.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  III  at  15  Hz.  Source  depth  100  m,  hydrophone 
depth  50  m.  Ice  roughness  3  m  RMS. 
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CAPTIONS  FOR  FIGURES  (cont'd)  ' 


Figure  32.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  III  at  15  Hz.  Source  depth  150  m,  hydrophone 
depth  50  m.  Ice  roughness  3  m  RMS. 


Figure  33.  Average  propagation  los3  as  a  function  of  range  for 
model  of  Table  III  at  15  Hz.  Source  50  m,  hydrophone 
depth  50  m.  Ice  roughness  3  m  RMS.  Averaging  inter-' 
val  11.5  km. 

Figure  34.  Average  propagation  loss  as  a  function  of  range  for 

model  of  Table  III  at  15  Hz.  Source  100  m,  hydrophone 
depth  50  m.  Ice  roughness  3  m  RMS.  Averaging  inter¬ 
val  11.5  km. 

Figure  35.  Average  propagation  loss  as  a  function  of  range  for 

model  of  Table  III  at  15  Hz.  Source  150  m,  hydrophone 
depth  50  m.  Ice  roughness  3  m  RMS.  Averaging  inter¬ 
val  11.5  km. 

Figure  36.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  VI  at  15  Hz.  Source  depth  50  m,  hydrophone 
depth  50  m.  Ice  roughness  3  m  RMS. 


Figure  37.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  VI  at  30  Hz.  Source  depth  50  m,  hydrophone 
depth  50  m.  Ice  roughness  3  m  RMS. 


Figure  38.  Propagation  loss  as  a  function  of  range  for  model  of 
Table  III  at  15  Hz.  Source  depth  50  m,  hydrophone 
depth  50  m.  Ice  roughness  4  m  RMS. 


Figure  39.  Average  propagation  loss  as  a  function  of  range  for 
model  of  Table  III  at  15  Hz.  Source  depth  50  m,  hy¬ 
drophone  depth  50  m.  Ice  roughness  4  m  RMS.  Averag¬ 
ing  interval  11.5  km. 

Figure  40.  Average  propagation  loss  as  a  function  of  range  for 
model  of  Table  III  at  20  Hz.  Source  depth  50  m,  hy¬ 
drophone  depth  50  m.  Ice  rouchness  3  m  RMS.  Averag¬ 
ing  interval  8.6  km. 

Figure  41.  Comparison  of  propag.-ion  loss  of  field  data  with  com 
puted  curve  of  Figure  40.  Predictions  by  ray  theory 
also  shown. 


Figure  A-l.  Layered  model 
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INTRODUCTION 

The  two  feature*  peculiar  to  the  polar  environment  that  most 
strongly  influence  underwater  sound  are  the  permanent  ice  cover 
and  the  velocity  structure  in  the  water.  Ice  movement  generates 
ambient  noise  and  the  ice  modifies  propagation,  particularly,  at 
high  frequencies,  by  scattering  waves  from  the  rough  ice  boundaries. 
Sound  velocity  is  a  function  of  temperature,  salinity,  and  pressure. 
The  relationship  among  these  variables  in  the  central  Arctic  Ocean 
is  such  that  sound  velocity  is  generally  an  increasing  function  of 
depth  from  the  surface  to  the  bottom.  Such  a  velocity  profile  is 
found  only  in  polar  waters.  The  sound  velocity  structure  is  very 
uniform  both  as  a  function  of  location  and  time  of  year.  Sounds 
are  transmitted  to  great  ranges  in  this  natural  Arctic  waveguide  or 
sound  channel  by  upward  refraction  in  the  water  and  repeated  reflec¬ 
tion  from  the  ice  canopy.  Typical  ray  paths  are  shown  in  Fig.  1.  A 
two-pound  explosion  of  TNT  has  been  heard  at  rrnges  exceeding  1,000 
km  (700  miles) . 

The  surface  sound  channel  of  the  Arctic  is  the  polar  extension 
of  the  deep  sound  channel  or  SOFAR  channel  of  the  nonpolar  oceans, 
but  the  Arctic  signals  are  often  quite  different  from  those  observed 
in  the  deep  channel,  largely  because  of  the  predominance  of  low-fre¬ 
quency  waves  in  the  Arctic.  A  regular  oscillatory  wave  train  at 
long  ranges  from  an  underwater  explosion  is  the  result  of  interference 
of  sounds  traveling  along  the  various  paths  between  source  and  detec¬ 
tor  (Kutschale,  1961,  1969).  The  exact  solution  of  the  wave  equation 
for  propagation  in  multilayered  media  in  terms  of  normal  modes  provides 
a  useful  method  for  describing  these  signals  in  detail  (Kutschale, 
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1969,  1970).  Figure  2  shows  a  typical  signal  with  waves  corresponding 
to  the  first  and  second  normal  modes.  On  Figure  3  a  typical  sound 
spectrogram  is  compared  with  a  computed  one.  The  agreement  between 
the  observed  and  computed  group-velocity  dispersion  is  excellent.  Fig¬ 
ures  2  and  3  show  that  the  low-frequency  waves  traveling  in  the  sound 
channel  are  coherent  over  long  ranges.  Thus  a  coherent  summation  of 
waves  is  appropriate  for  computing  propagation  loss  from  low-frequen¬ 
cy  continuous  wave  (cw)  sources. 

A  rapid,  accurate  method  is  developed  in  Appendix  A  of  computing 
propagation  loss  as  a  function  of  range  in  the  ice  covered  Arctic  0- 
cean.  Important  input  parameters  to  the  propagation  model  are  wave 
frequency,  source  and  detector  depth,  ice  roughness,  bottom  topogra¬ 
phy,  and  the  velocity  structure  as  a  function  of  depth  in  the  ice, 
water,  and  bottom  sediments.  Computational  speed  is  of  the  utmost 
importance  to  evaluate  effects  of  variations  of  these  parameters  on 
propagation  los3  as  a  function  of  range  since  a  large  number  of  mo¬ 
dels  must  be  considered. 

Such  a  rapid,  accurate  computational  method  is  the  Fast  Field 
Program  (FFP)  technique  for  wave  theory  introduced  by  Marsh  (1967), 

The  FFP  was  implemented  for  computation  on  a  digital  computer  by  Di- 
Napoli  (1971)  for  propagation  in  an  all-liquid  waveguide.  In  the 
Arctic  Ocean  solid  layers  as  well  as  liquid  layers  must  be  consid¬ 
ered,  and  a  convenient  starting  point  for  the  FFP  is  the  integral 
solution  of  the  wave  equation  derived  by  matrix  methods  in  Appendix 
A  for  propagation  from  point  harmonic  sources  in  a  multilayered 
liquid-solid  half  space. 

in  the  FFP  technique  the  integral  solution  of  the  wave  equation 
is  evaluated  rapidly  as  a  function  of  range  by  numerical  integration 
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employing  the  Fast  Fourier  Transform  (FFT).  Singularities  in  the  in¬ 
tegrand  corresponding  to  the  normal-mode  poles  are  removed  from  the 
axis  of  integration  by  including  attenuation  coefficients  for  com- 
pressional  and  shear  waves  in  each  layer.  In  a  liquid  layer,  of 
course,  both  the  shear  velocity  and  corresponding  attenuation  co¬ 
efficient  are  zero. 

Two  computer  programs  were  written  in  Fortran  IV  to  evaluate  the 
liquid-and  solid-bottom  integral  solutions  derived  in  Appendix  A.  A 
comparison  of  the  FFP  computations  with  those  computations  by  normal¬ 
mode  theory  from  the  corresponding  integral  solution  are  identical, 
but  the  FFP  technique  is  far  more  convenient  and  at  least  an  order  of 
magnitude  faster,  since  the  computations  are  done  directly  from  the 
integral  solution  without  first  computing  the  roots  of  the  dispersion 
equation  and  then  summing  the  normal  modes. 

In  the  multilayered  computations,  the  variation  of  compression- 
al  and  shear  velocities  and  density  with  depth  in  the  ice,  ocean,  and 
sediments  is  represented  by  a  series  of  plane  parallel  layers,  each 
with  constant  velocities  and  density  over  the  layer  thickness.  The 
last  layer  of  the  laminated  half  space  is  of  infinite  thickness.  The 
number  of  layers  is  adjusted  until  no  further  significant  change  of 
computed  results  is  obtained  by  a  finer  partitioning  of  the  observed 
velocity  and  density  data.  In  practice,  20  to  30  layers  are  usually 
sufficient  to  represent  the  velocity  and  density  data  in  the  deep  o- 
cean  over  the  frequency  ban'9  from  10  to  60  Hz.  At  higher  frequencies 
more  layers  may  be  necessary. 

Comparison  with  Normal-Mode  Theory 

The  FFP  approach  integrates  numerically  the  integral  solution  of 
the  wave  equation  derived  from  a  laminated  half  space.  For  guided 
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real  wave  number,  k,  axle  corresponding  to  the  normal-mode  poles. 

If  attenuation  coefficients  are  introduced  in  each  layer,  the 
poles  are  removed  from  the  real  k  axis  and  the  integration  can  be 
done  numerically.  The  effect  on  propagation  loss  as  a  function  of 
range  of  attenuation  in  each  layer  can  be  removed  after  the  integral 
is  computed  so  that  the  iinal  result  is  the  same  as  that  computed  by 
normal-mode  theory  with  no  attenuation  present  in  each  layer. 

Table  I  lists  the  parameters  of  a  model  from  which  computations 
were  done  by  normal-mode  theory  and  the  FFP  (8192  points  for  the  FFT) . 
The  velocity  profile  is  shown  in  Figure  4.  This  shallow-water  model 
exhibits  the  principal  features  of  normal-mode  propagation  at  low 
frequencies.  It  was  chosen  to  speed  the  numerical  work,  since  only 
between  two  and  five  normal  modes  propagate  in  the  band  from  10  to  25 
Hz.  Figures  5  to  8  show  the  comparisons  by  both  methods  of  propaga¬ 
tion  loss  as  a  function  of  range  at  10  and  25  Hz.  For  the  computa¬ 
tions  by  the  FFP  the  effect  on  propagation  loss  of  the  attenuation  in 
each  layer  was  removed  after  the  integral  wos  computed  so  that  the 
computations  from  the  integral  solution  by  both  methods  should  be  i- 
dentical.  That  this  is  the  case  is  apparent  from  Figures  5  to  8  at 
10  and  25  Hz. 

At  10  Hz  two  normal  modes  propagate  in  the  layered  system  and  at 
25  Hz  five  normal  modes.  The  absolute  value  of  the  integrand  at  10 
and  25  Hz  is  shown  in  Figures  9  and  10.  The  peaks  in  the  integrands 
correspond  to  the  wave  numbers  computed  by  normal-mode  theory.  The 
first  mode  corresponds  to  the  largest  k.  Progressively  higher  modes 
correspond  to  progressively  lower  k  values. 

Computation  of  Propagation  Loss 

This  section  presents  computations  of  propagation  loss  as  a  func¬ 
tion  of  ranrre  in  the  liaan  water  of  the  central  Arctic  fVacn 
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All  integrations  were  made  by  the  FFT  with  8192  points.  Great  care 
was  exercised  to  prevent  aliasing  in  the  FFT  computations. 

The  following  effects  of  the  waveguide  on  propagation  loss  are 
investigated* 

1)  An  ice  layer  at  the  surface 

2)  Rigidity  in  the  bottom  sediments 

3)  Source  depth 

4)  Small  changes  in  the  sound-velocity  profile 
in  the  upper  400  m  of  water 

5)  Variations  of  ice  roughness  along  the  path. 

Table  II  gives  a  layered  model  for  the  deep  Arctic  Ocean.  This 
model  represents  an  average  layering  in  the  ocean  based  on  about 
twenty  individual  sound  velocity  profiles  observed  at  various  lo¬ 
cations  and  at  different  times  of  the  year.  The  ocean  velocity  pro¬ 
file  is  partitioned  into  23  layers.  A  finer  partitioning  would  have 
a  negligible  effect  on  the  computations  over  the  frequency  range  frofn 
10  to  60  Hz.  The  4  km  water  depth  represents  propagation  over  the 
deep  abyssal  plains,  such  as  the  Canada  Abyssal  Plain. 

Figures  11  to  15  show  propagation  loss  as  a  tunction  of  range 
computed  rrom  tne  model  of  Tame  II.  The  fluctuating  propagation  loss 
from  tne  cw  source  in  the  band  from  10  to  60  Hz  is  in  response  to  the 
multipath  interference  of  wavos  traveling  between  source  and  detector. 
The  Integra'  ions  were  carried  out  over  a  range  of  wave  numbers  to  in¬ 
due  j  all  the  normal  modes  as  well  as  waves  from  which  energy  is  re¬ 
fracted  into  the  bottom  (leaky  modes).  The  contribution  of  the  leaky 
modes  to  propagation  loss  is  seen  as  a  high-frequency  oscillation  at 
short  ranges.  At  long  ranges  these  waves  have  been  eliminated  because 
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of  loss  of  energy  to  the  bottom.  The  effect  of  attenuation  in  each 
layer  was  removed  from  the  computations  of  Figures  11  to  15,  as  well 
as  in  all  subsequent  computations,  since  absorption  of  sound  in  the  wa¬ 
ter  is  negligible  at  the  very  low  frequencies  of  the  computations  over 
the  range  intervals  used.  However,  a  range  dependent  attenuation  for 
a  rough  ice  surface  is  included  by  multiplying  the  range  dependence 
of  pressure  by  the  formula  of  Mellen  and  Marsh  (1965). 

In  Figures  16  to  20  the  computations  are  repeated  for  the  layer¬ 
ing  in  the  water  of  Table  II,  but  a  solid  ice  sheet  3  m  thick  is 
placed  at  the  surface.  The  model  is  given  in  Table  III.  The  effect 
of  this  solid  ice  sheet  on  propagation  loss  as  a  function  of  range  ap¬ 
pears  to  be  negligible  over  the  frequency  range  10  to  60  Hz. 

Figures  21  to  30  show  the  effect  of  a  solid  bottom  on  propagation 

loss.  The  models  are  given  in  Tables  IV  and  V.  The  layering  in  the 

water  is  the  same  as  the  model  of  Table  II,  but  a  shear  velocity  of 
200  m/sec  is  included  in  the  bottom  sediments.  Such  low  shear  veloci¬ 
ties  are  commonly  measured  in  ocean-bottom  sediments  of  deep  basins 
(Sykes  and  Oliver,  1964).  It  is  apparent  that  the  effect  on  propaga¬ 
tion  loss  of  a  solid  bottom  is  small  if  the  ice  is  neglected.  When 
the  ice  sheet  is  included  at  the  surface  there  is  apparently  wave 
coupling  between  the  solid  bottom  and  the  solid  surface  layer  observed 
as  a  high-frequency  oscillation  of  propagation  loss  with  range.  How¬ 
ever,  the  average  propagation  loss  appears  to  be  nearly  the  same  as 
in  the  previous  cases. 

Figures  31  to  35  show  the  effect  of  source  depth  on  propagation 

at  15  Hz.  The  effect  of  source  depth  on  propagation  loss  is  more  ap¬ 

parent  if  the  propagation  loss  for  a  cw  source  is  smoothed  with  range. 
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This  smoothing  was  done  by  taking  a  running  average  of  propagation 
loss  in  range  computed  at  the  midpoint  for  each  range  interval.  This 
was  done  for  convenience  rather  than  taking'  the  running  average  with 
range  of  intensity  and  converting  this  to  propagation  loss.  The  lat- 

t 

ter  approach  is  rigorously  correct,  but  test  cases  showed  that  at 
15  and  20  Hz  the  difference  between  the  two  methods  of  averaging  was 
generally  less  than  three  dB.  There  is  only  a  slight  decrease  of  loss 
with  range  for  the  deeper  sources. 

Figures  36  and  37  show  the  effect  on  propagation  loss  at  15  and 
30  Hz  of  small  changee  in  stratification  in  the  upper  400  ro  of  water. 
The  model  is  given  in  Table  VI,  and  it  corresponds  to  a  single  meas¬ 
ured  profile.  It  is  apparent  by  comparing  Figures  36  and  37  with 
Figures  17  and  19  that  the  detailed  velocity  structure  in  the  upper 
layers  has  a  significant  effect  on  t^e  fluctuations  of  the  propaga¬ 
tion  loss  with  range,  but  the  average  loss  as  a  function  of  range  ap¬ 
pears  to  be  about  the  same  in  both  cases. 

Figures  38  and  39  show  the  effect  on  propagation  loss  of  a  rough¬ 
er  ice  sheet  than  used  in  the  previous  computations.  The  effect  on 
propagation  loss  of  this  ruugher  ice  sheet  is  not  apparent  when  com¬ 
pared  with  Figures  17  and  33,  but  it  would  be  at  longer  ramjres  beyond 
500  km  or  at  higher  frequencies  over  these  same  range  intervals. 

The  effects  of  variable  and  rough  bottom  topography  on  propaga¬ 
tion  loss  as  a  function  of  range  were  investigated,  but  not  in  detail. 
The  method  of  approach  is  to  assume  that  waves  passing  over  all  bottom 
topography  without  striking  the  bottom  predominate  in  the  signal.  Such 
an  assumption  is  substantiated  by  experimental  data.  This  is  imple¬ 
mented  in  the  computations  by  two  methods,  in  the  first,  the  layered 
lection  is  formed  for  the  deepest  portion  of  the  propagation  path.  In 
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layers  lying  above  all  topography,  low  attenuation  coefficients  are 
used  in  each  layer.  For  the  deeper  layers  much  higher  attenuation  co¬ 
efficients  are  used.  The  effect  of  this  high  attenuation  is  to  sup¬ 
press  waves  with  phase  velocities  greater  than  the  speed  of  sound  at 
the  depth  of  the  shallowest  part  of  the  bathymetric  profile.  The  sec¬ 
ond  method  is  to  integrate  only  over  the  range  of  wave  numbers  corres¬ 
ponding  to  waves  passing  over  all  bottom  topography.  This  method  has 
the  advantage  of  reducing  computer  time,  since  the  integral  is  comput¬ 
ed  over  a  shorter  range  of  wave  numbers. 

From  the  illustrations  it  is  apparent  that  an  ice  sheet  at  the 
surface,  a  low  rigidity  bottom,  or  small  changes  in  the  velocity  pro¬ 
file  mainly  affect  fluctuations  of  propagation  loss  with  range,  but 
do  not  appear  to  alter  significantly  average  values  of  loss  as  might 
be  measured  with  explosive  charges.  As  long  as  one  is  interested  on¬ 
ly  in  average  propagation  loss,  the  computations  for  an  all-liquid 
guide  appear  to  be  satisfactory.  On  the  other  hand,  if  geophones  are 
used  as  listening  devices,  the  ice  sheet  must  be  included  as  a  solid 
layer  at  the  surface  since  the  waves  are  elliptically  polarized  in 
the  ice  in  the  plane  of  propagation.  Particle  motion  for  guided  hy¬ 
droacoustic  waves  is  retrograde  elliptical  at  the  surface  and  pro¬ 
grade  elliptical  at  the  bottom  of  the  ice  (Kutschale,  1972).  The  pres¬ 
ent  theory  has  been  extended  to  compute  propagation  loss  to  geophones 
on  the  ice  surface  or  geophones  buried  at  depth  in  the  ice  (Kutschale, 
in  preparation). 

Comparison  with  Field  Data 

In  this  section  field  data  of  propagation  loss  as  20  Hz  are  com¬ 
pared  with  the  computed  propagation  loss.  The  comparisons  are  made 
for  an  average  propagation  loss  as  a  function  or  range  for  both  ex¬ 
periment  and  theory.  MeZlen  and  Marsh  (1965)  and  Buck  (1968)  have 


Page  14 


compiled  average  experimental  curves.  The  measurements  by  Mellen  and 
Marsh  and  Buck  lump  together  various  source  depths  between  40  and 
200  m  and  propagation  paths.  This  probably  accounts  for  the  spread 

of  measurements  at  each  frequency. 

Since  the  measurements  were  made  with  sxplosive  charges  analyzed 
in  various  filter  band  widths  reduced  to  a  one  Hz  band,  it  is  diffi¬ 
cult  to  compare  them  with  computations  from  a  cw  source.  A  meaning¬ 
ful  basis  for  comparison  is  the  smoothed  propagation  loss  as  a  func¬ 
tion  of  range  computed  from  a  running  average,  such  an  average  curve 
at  20  Hz  is  shown  in  Figure  40  for  the  model  of  Table  III.  Since  the 
preceding  computations  show  that  the  average  loss  is  not  affected  rou^h 
by  a  low  shear  velocity  in  the  bottom,  it  was  neglected.  Field  obser¬ 
vations  indicate  that  a  root-mean-square  (rms)  ice  roughness  of  3  m 
is  appropriate  over  the  paths  of  the  experiments. 

In  Figure  41  the  average  measured  curves  of  Mellen  and  Marsh  and 
Buck  are  compared  with  the  average  computed  one.  The  agreement  be¬ 
tween  the  two  sets  of  curves  is  quite  good,  A  detailed  comparison  be¬ 
tween  experiment  and  theory  would  require  close  control  of  the  experi¬ 
mental  situation.  Such  experiments  were  made  by  Kutschale  over  the 
flat  bottom  of  the  Abyssal  Canada  Plain  at  ranges  up  to  1300  km  by 
keeping  the  source  and  detector  depths  constant.  These  data  will  be 
analyzed  for  propagation  loss,  and  they  will  be  useful  for  a  detailed 
comparison  with  the  present  theory. 

Comparison  with  Ray  Theory 

Ray  theory  is  commonly  used  to  predict  propagation  loss  in  the 
ocean,  and  it  was  used  by  Aerophysics  Research  Corporation  in  their 
work.  If  only  propagation  loss  at  a  single  range  point  on  a  path 
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is  required,  it  may  be  faster  to  use  ray  theory  than  the  direct  in¬ 
tegration  (FFP).  The  question  is  often  raised,  however,  about  the 
reliability  of  predictions  by  ray  theory  at  very  low  frequencies. 

To  attempt  to  answer  this  question,  the  FFP  computations  are  compared 
with  the  ray  predictions  in  Figure  41.  Propagation  loss  by  ray  the¬ 
ory  was  computed  by  summing  the  intensities  for  each  ray  incoherent¬ 
ly  (phase  was  neglected).  These  computations  should  thus  be  compared 
with  the  smoothed  propagation  loss  computations  by  the  FFP.  it  is  seen 
from  Figure  41  that  at  ranges  beyond  100  km  the  predictions  by  both 
methods  agree  reasonably  well.  At  ranges  less  than  100  km  the  two 
sets  of  computations  diverge.  The  ray  theory  predictions  of  loss  are 
too  low  compared  both  to  the  field  data  and  the  FFP  results.  The  rea¬ 
sonably  good  agree. nent  at  long  ranges,  however,  is  gratifying  and  sup¬ 
ports  the  application  of  ray  theory  for  the  intensity  computations  at 
Aerophysics  Research  Corporation  at  ranges  exceeding  50u  km. 
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Layer 

TABLE  I 

Layered  Model,  water  depth  350  m 

Compreseional 

velocity, 

m/sec 

Shear 

velocity, 

m/sec 

Density,  gm/cro 

1 

1434.0 

0.0 

1.03 

2 

1437.0 

0.0 

1.03 

3 

1440.0 

0.0 

1.03 

4 

1443.0 

0.0 

1.03 

5 

1446.4 

o.o 

1.03 

6 

1450.0 

0.0 

1.03 

7 

1453.2 

0.0 

1.03 

8 

1600.0 

0.0 

1.03 

Layer 

thickneee, 

m 


50.0 

50.0 

50.0 

50.0 

50.0 

50.0 

50.0 

OO 
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TABLE  II 

Layered  Model,  water  depth  4000  ro.  liquid  bottom 


Layer 

Conpreeslonal 

velocity, 

m/sec 

Shear 

velocity, 

w/sec 

3 

Density,  gm/cm 

Layer  thickness 
in 

1 

1434.0 

0.0 

1.03 

50.0 

2 

1437.0 

0.0 

1.03 

50.0 

3 

1440.0 

0.0 

1.03 

50.0 

4 

1443.0 

0.0 

1.03 

50.0 

5 

1446.4 

0.0 

1.03 

50.0 

6 

1450.0 

0.0 

1.03 

50.0 

7 

1453.2 

0.0 

1.03 

50.0 

8 

1456.0 

0.0 

1.03 

250.0 

9 

1460.0 

0.0 

1.03 

200.0 

10 

1463.0 

0.0 

1.03 

200.0 

11 

1467.0 

0.0 

1.03 

280.0 

12 

1472.0 

0.0 

1.03 

220.0 

13 

1476.0 

0.0 

1.03 

260.0 

14 

1480.0 

0.0 

1.03 

240.0 

15 

1484.0 

0.0 

1.03 

280.0 

16 

1489.0 

0.0 

1.03 

320.0 

17 

1494.0 

0.0 

1.03 

240.0 

18 

1498.0 

0.0 

1.03 

200.0 

19 

1502.0 

0.0 

1.03 

200.0 

20 

1506.0 

0.0 

1.03 

280.0 

21 

1511.0 

0.0 

1.03 

240.0 

22 

1515.0 

0.0 

1.03 

240.0 

23 

1600.0 

0.0 

1.70 
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TABLE  III 


Layer 

Levered  Model. 

water  depth  plus 

Ice  thickness 

4000  m. 

Compressional 
velocity,  n/aec 

licruid  bottom 

Shear 

velocity,  , 

m/sec  Density,  gin.  cnr 

Layer 

thickness,  m 

1 

3500.0 

1800.0 

0.90 

3.0 

2 

1434.0 

0.0 

1.03 

47.0 

3 

1437.0 

0.0 

1.03 

50.0 

4 

1440.0 

0.0 

1.03 

50.0 

5 

1443.0 

0.0 

1.03 

50.0 

6 

1446.4 

0.0 

1.03 

50.0 

7 

1450.0 

0.0 

1.03 

50.0 

8 

1453.2 

0.0 

1.03 

50.0 

9 

1456.0 

0.0 

1.03 

250.0 

10 

1460.0 

o.o 

1.03 

200.0 

11 

1463.0 

0.0 

1.03 

200.0 

12 

1467.0 

0.0 

1.03 

280.0 

13 

1472.0 

0.0 

1.03 

220.0 

14 

1476.0 

0.0 

1.03 

260.0 

15 

1480.0 

0.0 

1.03 

240.0 

16 

1484.0 

0.0 

1.03 

280.0 

17 

1489.0 

0.0 

1.03 

320.0 

18 

1494.0 

0.0 

1.03 

240.0 

19 

1498.0 

0.0 

1.03 

200.0 

20 

1502,0 

0.0 

1.03 

200.0 

21 

1506.0 

0.0 

1.03 

280.0 

22 

1511.0 

0.0 

1.03 

240.0 

23 

1515.0 

0.0 

1.03 

240.0 

24 

1600.0 

0.0 

1.70 
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TABLE  IV 


Layer 

Layered  Model,  water  depth  4000  n..  solid 

bottom 

Compreseional 

velocity, 

m/sec 

Shear 

velocity, 

m/sec 

Den»ity,gm/ctn^ 

Layer  thickness 
m 

1 

1434.0 

0.0 

1.03 

50.0 

2 

1437.0 

0.0 

1.03 

50.0 

3 

1440.0 

0.0 

1.03 

50.0 

4 

1443.0 

0.0 

1.03 

50.0 

5 

1446.4 

0.0 

1.03 

50.0 

6 

1450.0 

0.0 

1.03 

50.0 

7 

1453.2 

0.0 

1.03 

50.0 

8 

1456.0 

0.0 

1.03 

250.0 

9 

1460.0 

0.0 

1.03 

200.0 

10 

1463.0 

0.0 

1.03 

200.0 

11 

1467.0 

0.0 

1.03 

280.0 

12 

1472.0 

0.0 

1.03 

220.0 

13 

1476.0 

0.0 

1.03 

260.0 

14 

1480.0 

0.0 

1.03 

240.0 

15 

1484.0 

0.0 

1.03 

280.0 

16 

1489.0 

0.0 

1.03 

320.0 

17 

1494.0 

0.0 

1.03 

240.0 

18 

1498.0 

0.0 

1.03 

200.0 

19 

1502.0 

0.0 

1.03 

200.0 

20 

1506.0 

0.0 

1.03 

280.0 

21 

1511.0 

0.0 

1.03 

240.0 

22 

1515.0 

0.0 

1.03 

240.0 

23 

1600.0 

200.0 

1.70 
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TABLE  V 

Laytrad  Modal,  water  plua  lea  dapth  4000  m. 


aolld  bottom 


Layer 

Compreaaional 

velocity, 

tn/eec 

Shear 

velocity, 

m/sec 

3 

Density,  gtn/cro 

Layer  thickness 
m 

1 

3500.0 

1800.0 

0.90 

3.0 

2 

1434.0 

0.0 

1.03 

47.0 

3 

1437.0 

0.0 

1.03 

50.0 

4 

1440.0 

0.0 

1.03 

50.0 

5 

1443.0 

0.0 

1.03 

50.0 

6 

1446.4 

0.0 

1.03 

50.0 

7 

1450.0 

0.0 

1.03 

50.0 

8 

1453.2 

0.0 

1.03 

50.0 

9 

1455.0 

o.o 

1.03 

250.0 

10 

1460.0 

0.0 

1.03 

200.0 

11 

1463.0 

0.0 

1.03 

20C.0 

12 

1467.0 

0.0 

1.03 

280.0 

13 

1472.0 

0.0 

1.03 

220.0 

14 

1476.0 

0.0 

1.03 

260.0 

15 

1480.0 

0.0 

1.03 

240.0 

16 

1484.0 

0.0 

1.03 

280.0 

17 

1489.0 

0.0 

1.03 

320.0 

18 

1494.0 

0.0 

1.03 

240.0 

19 

1498.0 

0.0 

1.03 

200.0 

20 

1502.0 

0.0 

1.03 

200.0 

21 

1506.0 

0.0 

1.03 

280.0 

22 

t 

1511.0 

0.0 

1.03 

240.0 

23 

1515.0 

0.0 

1.03 

240.0 

24 

1600.0 

200.0 

1.70 

© o 
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TABLE  VI 


Layer 

Layered  Model,  water  plus  ice  depth  4000  in 

# 

,  Layer 

thic’cne 
m 

Coinpressional 
velocity,  m/sec 

licruid  bottom 
Shear  velocity, 
m/sec 

Density, gm/cm 

1 

3500.0 

1800.0 

0.90 

3.0 

2 

1435.2 

0.0 

1.03 

37.0 

3 

1438.0 

0.0 

1.03 

10.0 

4 

1438.0 

0.0 

1.03 

30.0 

5 

1441.0 

0.0 

1.03 

60.0 

6 

1443.0 

0.0 

1.03 

40.0 

7 

1445.0 

0.0 

1.03 

40.0 

8 

1449.0 

0.0 

1.03 

40.0 

9 

1453.0 

0.0 

1.03 

60.0 

10 

1455.0 

0.0 

1.03 

40.0 

11 

1457.0 

0.0 

1.03 

100.0 

12 

1453.5 

0.0 

1.03 

140.0 

13 

1460.0 

0.0 

1.03 

200.0 

14 

1463.0 

0.0 

1.03 

200.0 

15 

1467.0 

0.0 

1.03 

280.0 

16 

1472.0 

0.0 

1.03 

220.0 

17 

1476.0 

0.0 

1.03 

260.0 

18 

1480.0 

0.0 

1.03 

240.0 

19 

1484.0 

0.0 

1.03 

280.0 

20 

1489.0 

0.0 

1.03 

320.0 

21 

1494.0 

0.0 

1.03 

240.0 

22 

1498.0 

0.0 

1.03 

200.0 

23 

1502.0 

0.0 

1.03 

200.0 

24 

1506.0 

0.0 

1.03 

280.0 

25 

1511.0 

0.0 

1.03 

240.0 

26 

1515.0 

0.0 

1.03 

240.0 

27 

1600.0 

0.0 

1.70 

oo 

+769.000  SEC 
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Fig.  11 
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DEFINITION  OF  SYMBOLS  FOR 
APPENDIX  A  AND  APPENDIX  B 

cotnpressional-wava  velocity  in  the  m-th  layer 
shear-wave  velocity  in  the  m-th  layer 

compressiona 1-wave  attenuation  coefficient  in  the  m-th  layer 

shear  wave  attenuation  coefficient  in  the  m-th  layer 

thickness  of  the  m-th  layer 

vertical  coordinate 

range  between  source  and  detector 

time 

angular  frequency 
phase  velocity 
wave  number 

r? 

real  part  of  a  complex  number 

imaginary  part  of  a  complex  number 

magnitude  of  a  complex  number  £ 

velocity  potential  in  the  m-th  layer 

normal  stress  in  the  m-th  layer  parallel  to  z  axis 
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t  tangential  streaa  in  the  m-th  layer 

JD  t  pressure  In  the  m-th  liquid  layer 
'AW. 

*  horizontal  particle  displacement  in  the  p-th  layer  in  the 
radial  direction 


vertical  particle  displacement  in  the  m-th  layer 
horizortal  particle  velocity  in  the  m-th  layer  in  the  radial 


direction 

t  vetical  particle  velocity  in  the  m-th  layer 


density  in  the  m-th  layer 
density  at  the  source 
Bessel  function  of  order  0 

Y  Bessel  function  of  order  0 

Hankel  function  cf  the  second  kind  of  order  0 
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APPENDIX  A 

INTEGRAL  SOLUTION  FROM  POINT  HARMONIC  S0URCE9 

The  solution  of  the  wave  equation  presented  here,  based  on 
the  Thomson-Haskell  matrix  method  (Thomson,  1950?  Haskell,  1953), 
follows  Harkrider  (1964)  for  the  solution  of  the  wave  equation 
in  an  n-layered  solid  half  space.  Layer  matrices  of  the  type 
given  by  Dorman  (1962)  for  computing  dispersion  in  an  n-layered 
liquid-solid  half-space  are  used  for  the  liquid  layers.  An  ap¬ 
plication  of  the  theorem  that  the  inverse  of  the  product  matrix 
for  the  layered  system  above  the  source  has  the  same  form  as  the 
inverse  of  a  layer  matrix  reduces  the  integrand  of  the  integral 
solution  to  a  simple  form  in  terms  of  elements  of  product  ma¬ 
trices  derived  in  the  source-free,  plane-wave  case. 

Consider  an  interbedded  stack  of  plane  parallel  liquid  and 
solid  layers  resting  on  a  solid  half  space.  The  layers  are  shown 


in  Figure  A-l. 
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Solutions  Jn  the  m-th  liquid  layer  of  vertical  particle  veloc¬ 
ity  and  stress  are 

C  J0(kr)e.tWtdk 
0  (1) 


We  divide  the  liquid  source  layer  and  detector  layer  each  in¬ 
to  two  layers  as  shown  in  Figure  A-l.  At  the  bottom  of  the  Sx  lay¬ 
er  the  pressure  is  continuous.  The  vertical  particle  velocity  is 
continuous  everywhere  in  the  plane  between  the  Sl  and  S2  layers  ex¬ 
cept  at  the  point  source  where  the  liquid  above  and  below  the  source 
moves  in  opposite  directions.  This  may  be  expressed  by  writing 


« 

^(^0  ~  2  K> 


In  the  liquid  source  layer  for  the  -dependent  particle  velo 


cities  and  stresses 


i 

A 

where  *■  the  horizontal  particle  velocity  at  the  bot 


tom  of  the  first  liquid  layer  above  the  half  space.  For  the  layers 
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and  for  the  layer*  above  the  s~  layer 


AV  °v  bV> . k 


The  matrices,  am,  for  liquid  layers  and  ,bm/  for  solid  layers  be¬ 
tween  liquid  layers  are  in  terms  of  4  x  4  matrices: 

l  o  0  O 

(X  -  0  (CW)ax  0 

o  CA-'wv')^  o 

oo  o  l 
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in  term*  of  element*  of  the  eolid  layer  matrix 


(A<wv\3  (cx**)^ 

(A*^9a 

CaAa 


The  matrix  element,  for  liquid  and  .olid  layer,  are  given  in 


Appendix  B. 


We  now  write  for  the  .olid  half-.pace  (Haskell,  1953) 


1 

r  A 

^  VV^| 

A  1  1  « 

^  W\  ~  Ay^ 

'lU 

\l 

• 

/\ 

<-< 

[ 

i  * 

J 

am,  •— 
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equal  zero  and  hence 


i 


A  _ 

1 

t 

-A. 

V 

*  E." 

A 

A 

<< 

> 

1 

(6) 
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We  may  prove  following  Harkrider  (1964)  that  the  inverse  of  the 
product  matrix  Ac  hao  the  same  form  as  the  inverse  of  a 
layer  matrix  and  hence 
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(11) 
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I* 


we 

■ 

v\* 

I 


define  J"  ~  E  A 


U) 


< 


”*1 


Eliminating  ^ 


and  write  Equation  (7) 


<\ 


4 


^  ^33  ^ 


~^4i  ^4-i  X 


'43  '*4'i 


*> 

'V/' 

X 

'f 

.2 

(  V  Tai^V/  +(Xa“  Xu.)X  +  (12) 

(X3  ~  Xy)Y  +(X^  ~  -  0 


and  eliminating 

(^1  ~^4iW  v  (^3a.-W)X  +  (i3) 

UM-V>'(  +  ( TM  -  T^Z  -  0 
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If  we  let 


*-•  *  3"^  Tai 


Cr  ~  Tj-j  -  Taj 


K- 

3h  -■ 

> 

N  - 

^4i 

M  - 

^■ja.  “ 

^4  a. 

M  « 

^*33  ~ 

T^<i 

S  - 

^34  “ 

■  T44 

and  solve 

Equations  (12)  and  (13) 

for  W  and  X  respectively, 

v/-  - 

K  Y  _ 

ux 

■fcV 

-B.Z- 

u  27 

X- 

-N.W 

_  H' 

1-^H- 

M 

M 

M 

(14) 

(15) 
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or  substituting  Equation  14  in  Equation  15 


and  solving  for  X 

Y=  (G-N-HL^  +  (RN-SQZ- 

ML-  NK 


and  therefore  from  Equations  (10)  and  (11) 

X  ~  t  & N -  H 

NK-ML 

or  from  the  Equation  9 


0 


(&N-UL)(A^3a 

NK-ML. 
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Hence,  the  integral  solution  for  UJ*  is  ,frora  Equation  1 

iiot  0 

dropping  the  0  term, 


/°°  (GN-  ^L)(A$,\a 
<  _ 

0  (NK-ML) 


Jo(kr)<3,  Kd k. 

If  a  hydrophone  or  vertical  particle  velocity  detector  is 


located  in  a  liquid  layer  at  the  bottom  of  the  layer  in  Fig¬ 
ure  A-l  . 


where 
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or 


(+»V'  •  i ' 


Hence,  the  integral  solution  for  pressure  at  the  bottom  of  layer 


D,  is 


•fe 


"C 


|(&N  -  HU)(AS) 


(nk-ml) 


I  ■ 


(16) 


(A4l.(!<r)Akclk 

which  is  the  desired  integral  solution  (  Equation  3. 

For  a  constant  pressure  source  of  amplitude  P  dynes/cm2  at 


1  m  we  multiply  Equation  16  by  0 

(.'•of* 


and  get 


1  c  Vps 


V0  I (&N-MuyA^) 

(nk-ml.) 


t-ml 


( ^ D,\a  kdk 
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If  the  bottom  is  liquid  the  corresponding  formula  is  (Kutschale 
1970)  in  terms  of  4  Joy  4  layer  matrices  of  Appendix  B 

D|  C  Wl 


A 


a  a  +  t  A^) 


(A^  <2  T0(J?r)  kdV; 


Carrying  out  the  matrix  multiplication 
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I 


and  eliminating  and 

L,  K.  G,  N,  M  and  H  are 


explicit  expressions  for 
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W«  now  introduce  attenuation  coefficient! 

) 

in  ®ach  layer  which  removes  the  singularities  in  the 

integrand  from  the  real  k  axis.  Thie  is  done  by  setting  |<  L 

*  \ 


complex  in  each  layer: 

r*  ^  i  d> 


k*\ 


Since 


c*o 

A  Ckr)-  J0(krW 

the  integral  solution  for  a  constant  power  source  may  be  written 

^“[(gn-hlXaA 

- —  -K\ 


(NK-  ML) 


(AXa  ^  ^  (  kr)  k  d  k 


if  we  remember  that  only  the  part  of  U^//i  V 

J.CHO 


orresponding  to 


is  desired. 
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Let 


and  set 


_  ^  )(G-  n-u  y 

^KNK-ML.) 


V. 


CAI>,)3a  5  ^ 


This  approximation  of 


is  valid  beyond  a  wave¬ 


length  from  the  source.  The  integral  solution  is  then 


-ckr 


dk. 


(17) 
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Let 


*«.  *  ro  v  n^*~  n^o,  -•■  L»-| 


with 


AkA r= ay 

4l 


must  be  an  integral  power  of  two. 
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Substituting 


in  Equation  17 


~L  7^  +^4kr„) 


Vi 


-atrc/wv-n. 


l/ni 


FCM  *  Re^FOO]  +- cT/*v[FCky\ 


or 


rfc. 


(Vyl-o 


This  formula  is  rapidly  evaluated  at  n  range  points  jfl 


the  FFT  for  the  complex  input  X  ■**v  Y 


n. 


by 


The  upper  limit  of  integration  is  chosen  to  include  waves  cor¬ 
responding  to  nonzero  values  of  the  integrand  in  the  layered  system. 
The  upper  limit  of  integration  is 


k0+  CN-l)Ak 
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from  which 

Ak  »  ^  ~  *° 

M  - 1 


and 


snr 

NAk 


Propagation  loss  at  range  ^”r\_  is  defined 

■pi*tO*-aoHo«  'ki 

'Mo  “p 
1  o 

If  the  attenuation  coefficients  are  equal  in  each  layer  the  effect 


on  propagation  loss  of  the  attenuation  is  removed  by  multiplying 
by 

l 


h 


W<l>  Wll«  W  V  IIMW  W  A  Wil  A  U  ~  z  *•  »  ^  »  •  ••3 

-Rro  -  e ' ^  ea/VJV" . 


In  this  case 


puo-.aodc* 


* 


\0 


T. 
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The  effect  of  attenuation  by  a  rough  ice  sheet  is  introduced 


in  the  propagation  loss  computations  by  multiplying 
by 


where 


4 

H 

Co 


wave  frequency 

RMS  ice  roughness  in  m 

surface  sound  velocity  in  the  water  in 
m/sec 

range  for  one  cycle  of  a  ray  of  phase 
velocity  c. 


This  formula  is  derived  by  Marsh  et  al  (1961)  and  was  used  by 
Mellen  and  Marsh  (1965)  for  their  work  on  Arctic  propagation,  it 
has  proved  useful  at  low  frequencies  if  Q  ,  and 
are  computed  from  a  ray  vertexing  at  1000  m  for  deep  water  propaga¬ 


tion. 
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APPENDIX  B 

LAYER  MATRIX  ELEMENTS 


Solid  Layers: 
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Liquid  layers: 

=  &»*\x  ~  0 &/w\  * 

-  O 

(Q^^-  L  ^  ^  1 

1/0  P*v 

lA*^  -  L  °°P-  ^>vv^ 
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